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Vortices in a thin-film superconductor interact logarithmically out to a distance on the order of the 
two-dimensional (2D) magnetic penetration depth Ax, at which point the interaction approaches a 
constant. Thus, because of the finite Ax, the system exhibits what amounts to an intrinsic finite 
size effect. It is not described by the 2D Coulomb gas but rather by the 2D Yukawa gas (2DYG). 
To study the critical behavior of the 2DYG, we map the 2DYG to the massive sine-Gordon model 
and then perform a renormalization group study to derive the recursion relations and to verify that 
Ax is a relevant parameter. We solve the recursion relations to study important physical quantities 
for this system including the renormalized stiffness constant and the correlation length. We also 
address the effect of current on this system to explain why finite size effects are not more prevalent 
in experiments given that the 2D magnetic penetration depth is a relevant parameter. 



I. INTRODUCTION 

Interest in the critical behavior of vortices in two- 
dimensional (2D) superconductors can be traced back to 
at least the studies of Kosterlitz and Thouless |j],E| and 
Berezinskii || in the early 1970's. It was stated in Rcf. |l]] 
however, that the theory as formulated did not apply 
to 2D superconductors where there is a coupling to the 
magnetic field. Because of this coupling, the interaction 
is logarithmic only up to the length scale Ax = 2A 2 /d, 
(where A is the London penetration depth and d is the 
film thickness,) at which point it approaches a constant. 
The application of the Kosterlitz-Thouless-Berezinskii 
(KTB) transition to 2D superconducting systems dates 
only to the late 1970's and the work of Beasley, Mooij, 
and Orlando Q and of Halperin and Nelson jq] who 
pointed out that Ax can be quite large (relative to the 
film size) in superconducting thin films and that there- 
fore these films should behave as the charged analog of 
superfluid thin films. As a result, a flurry of experimen- 
tal papers reported investigations of KTB behavior in 
thin-film superconductors in the early 1980's. 

Sparked partially by the discovery of the high- 
temperature superconductors, interest in the the critical 
behavior of 2D superconductors (SC's) and Josephson 
Junction arrays (JJA's) has been renewed in this decade. 
Sujani, Chattopadhyay, and Shenoy j^] have derived re- 
cursion relations for 2D superconductors that incorporate 
nonzero current drives. Minnhagen and cowor kers @ 
have been investigating the dynamics of the transition 
numerically while Ammirata et al. [ll] , |l2[ ] have applied 
dynamic scaling to I-V data from 2D superconductors 
and JJA's. 

The limiting large scale Ax , a length intrinsic to 2D su- 
perconductors, acts as an unavoidable source of finite size 
effects. Other sources, characterized by other lengths, 



are extrinsic, arising from the finite size (or extent) of 
a system and from a current. Some of the recent work 
on 2D SC's and JJA's has focused on the expected finite 
size effects (of any origin) in these systems. Simkin and 
Kosterlitz [ fl3| have investigated numerically and analyt- 
ically the finite size effects in these systems, while Repaci 
et al. jl4| claimed experimental evidence of finite size ef- 
fects in ultra thin film YBCO. Herbert et al. (l5) have 
extended this work to JJA's. We will focus here on the 
influence of a finite Ax. , which we refer to as the intrinsic 
finite size effect. 

Despite all of the work on 2D SC's and JJA's, there 
is still no analytical theory of the influence of intrinsic 
finite size effect on their critical properties. This is par- 
ticularly striking since it is generally accepted that the 
intrinsic finite size effect should be relevant, yet reports 
of its experimental observation are relatively sparse. Our 
objective in this paper is to fill in the theoretical gap and 
to address the small number of experimental reports of 
finite size effects. We perform a RG study of the vortices 
in 2D superconductors by modeling them as a 2D Yukawa 
gas (2DYG) and then mapping this system to the mas- 
sive sine-Gordon (SG) Hamiltonian. We then perform a 
momentum-space shell RG study to derive the recursion 
relations, which we then analyze to study the shift in the 
critical temperature and the behavior of the system near 
it. We will also address the effect of the presence of a 
driving current on the critical behavior of the system. 

This paper is organized as follows. In Section ||, the 
partition function for the 2DYG is given and shown to 
map to the massive sine-Gordon. In Section [II Oj , the 
recursion relations are derived. These equations are ana- 
lyzed in Section |l| and their experimental ramifications 
are explored in Section IV . Finally, the paper is summa- 
rized in Section 
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II. MODEL AND RENORMALIZATION GROUP 
STUDY 

As indicated above, our approach is to model the vor- 
tices in two dimensional superconductors as a vortex gas 
interacting via the 2D Yukawa potential. Therefore, we 
begin this section by describing the partition function 
for the 2DYG. We then will map it to the massive sine- 
Gordon model, and perform a RG study on the latter. 



A. Model 

The partition function for the neutral vortex gas (iV+ 
vortices and JV_ = N + antivortices) in a 2D supercon- 
ductor is 
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where 2N is the total number of particles (N — N±), 
are the coordinates of the zth charge pi — ±p and y — 
exp(/3/i) is the fugacity, where [i = —E c and j3 = X/ksT. 

V(R) is the interaction of a pair of vortices, the 
strength of which is p 2 = im 2D Ti 2 /2m (n 2D = n s d is 
the areal superfluid density, n s is the superfluid elec- 
tron density, and m is the mass of a free electron.) 
For the 2D Coulomb gas, the interaction is logarithmic 
V(R) = ln(R/r) (where r is the vortex core size) for all 
lengths. But for a superconductor, the interaction is only 
logarithmic out to Ax, at which point it approaches 
a constant as 1/R due to the coupling to the magnetic 
field. We will approximate this behavior for the potential 
by writing it in terms of a modified Bessel function, 



V(R) = ]m(X±/t) -K (R/Xa_), 



(2) 



which has the limits calculated by Pearl. [Q We will re- 
fer to the vortex gas in 2D superconductors as the 2D 
Yukawa Gas (2DYG) since this potential is the 2D equiv- 
alent of the more familiar 3D Yukawa potential. 



B. Mapping to the Massive sine-Gordon 

A standard way of treating the 2DCG is to map it 
to the 2D sine-Gordon Hamiltonian. |L6j For the case of 
the 2D vortex gas in a superconductor, as we now review, 
the correct mapping is to the massive sine-Gordon model: 



where X = \/[2n(3p 2 ] is the coupling strength or "stiff- 
ness," <j> is the SG field and m is the mass, which we will 
relate to X± below. For m — 0, Eq. (||) reduces to the 
SG Hamiltonian. 

In Ref. [p| the mapping of a layered vortex gas with 
arbitrary interactions to a SG-like Hamiltonian was stud- 
ied. By examining the 2D limit of this mapping, we 
can establish the link between the 2DYG and the mas- 
sive sine- Gordon. Following Ref. (TJ], one can show that 
Eq. (|l|) can be written as 

J gflr) exp { -H + 2y J d 2 r cos[</>(r)] } 

/ D</>(r) exp[-tf ] ' 1 > 

where the Hamiltonian H a is Gaussian in the field <f>. 
The connection between the sine-Gordon field and the 
potential V(R) is given by: 



I3 P 2 V(R) =< <f>(R)4>(0) > - < 2 (O) > , 



(5) 



where (...)o denotes an average with respect to a Hamil- 
tonian H Q Gaussian in the field <f>. 

Making use of Eq. (J|), one can readily show 
that V{R) = ln(A x /r) - K (R/X±) when H = 
1 J d 2 r{[V(j)(r)] 2 + [m0(r)] 2 }, (3p 2 = 1/[2ttX], and m = 
t/Ax- The calculation of the right hand side of Eq. (||) 
is straightforward and quite similar to that done for the 
layered case in Ref. (See the calculation of g s {R) 

in Appendix B of that reference.) After performing the 
Gaussian integral over the fields, one is left with 



Pp 2 v(R) = 



d 2 q 1 — cos(q • R) 
(2tt) 2 21{q 2 + m?) ' 
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where the angle is integrated from to 2ir and q is in- 
tegrated from zero out to its infrared cutoff qo ~ 2tt/t. 
Performing the integral over the angle first and then over 
q, one obtains 



PP 2 V(R) = _L{[1 - h(qoR)] InR 

-K (mR) -ln(mR)\, 



(7) 



where R is now expressed in units of r and q in units 
of 1/r. Neglecting the oscillating behavior (2(| of the 
zeroeth order Bessel function 3 (qoR) at large R, which 
arises from the sharp cutoff, this has the desired behavior 
ofEq. (|). 

This completes that mapping between the 2DYG and 
the massive 2D SG Hamiltonian. We will now proceed 
to a RG study on the latter. 
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C. Renormalization Group Study 

The advantage of the above mapping is that it is eas- 
ier to perform a momentum space RG study on Eq. (0) 
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than a real space RG calculation based on Eq. ([!]). In 
this subsection we show how to perform a RG study in 
momentum space on Eq. (^) to obtain the recursion re- 
lations for the three parameters in the system X, y, and 
to. The recursion relations for this system were previ- 
ously derived in Ref . H] , but as we will explain below, a 
technical point was neglected which rendered the results 
formally incorrect. 

To briefly overview the process, the RG study consists 
of integrating out the large wavevector q components in a 
thin shell: [qo(l — e) < q < qo], which correspond to small 
scale structure. The system is then shrunk back so that 
the cutoff in momentum space is restored to the original. 
We follow the procedure of Knops and den Ouden Jll| 
as outlined in Rcf. jl9| for a layered vortex gas. Because 
the calculation is closely related to the previous work of 
the present authors, |ll| we do not need to show much 
detail here. 

The first step in the RG study is to Fourier analyze 
the first term in the free energy, 



f = -\I ^m^-q)W 
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it is then convenient to define the Gaussian terms in 
Eq. (|) as 



F a = cf ) (q)^-q)[I(q 2 +m 2 )}. 
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and 
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Integrating Eqs. (14) and ([19) we have 
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and 



g(R) = -=. [(1 - J (q R))laR - K (mR) - ln(mR)} . 

(17) 

Substituting these into A(X, m) and expanding to first 
order in m, one finds 19 



A(I,m) = A(T) - A m (T)mlnm. 



(18) 



The values of A(T) and A m (T) depend upon the cutoff 
used and the nature of core used. Here, as in Ref. p9[ , 
we will assume a sharp cutoff in momentum space and a 
hard-core. 

The recursion relations can be put in more convenient 
differential form. Introducing also the stiffness constant, 
K = l/[27rl](= (3p 2 ), we can write the final result in the 
form 



dK 1 2 2/ , 
— — = — K y (1 — 4.5mmm), 
de 2 y K h 
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There are two differences between the 2DYG and 
the layered SG model of Ref. Q: First, F a for the 
layered case is not given by Eq. (^), but rather by 
4>k(q)4'-k{—q)[Zq 2 + 2K(1 — cos(fcs))] where s is the in- 
tcrlayer distance and k is the z-axis wave vector. (K also 
has a different meaning there.) Second, the system here 
is purely two-dimensional and so the integration over the 
third dimension k in Ref. (T^| does not appear. Keeping 
these two distinctions in mind and carrying out the RG 
study in the manner of Ref. |lj| , one finds 



K 



and 



where 
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A(I, m) = J d 2 R{q ■ H) 2 f(R)e- 9(R)/2 , (13) 
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and 



de 2 y 1 



dm , K 

— =m(l-y —). 
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Our definition of the parameter K differs from others in 
the literature by a factor of 2-k. 

The to = limit of these equations correctly reduce to 
the original Kosterlitz recursion relations for the 2DCG. 
Q Furthermore, the m corrections to those relations are 
of the form that one would expect physically. |22] For ex- 
ample, in the recursion relation for K, the screening ef- 
fect of small pairs on K is increased by m (since m « 1). 
This is because the vortex pairs are more weakly bound 
for a superconductor than for a superfluid and so they are 
more susceptible to screening. (This is quite consistent 
with the opposite result found in the case of a layered vor- 
tex gas. There, the vortex pairs are more strongly bound 
than in the 2DCG and the screening effects in the recur- 
sion relation for K diminished. p2] , p3[ |) In the recursion 
relation for y, the number of vortex pairs increases in the 
presence of finite size effects because the energy of a vor- 
tex pair is reduced. In the recursion relation for m, there 
are two terms. The first enters through the shrinking 
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scale step of the RG calculation. The second, of 0(y 2 ), 
enters through the integrating out of small scale struc- 
ture, which in this case consists of the small vortex pairs. 
This explains the y 2 dependence. This second term arises 
in a somewhat complicated way from the renormalization 
of the coefficient of the (V</>(r)) 2 term in Eq. @. This 
enters into the mass term m because the coefficient of 
(m<fi(r)) 2 is Im 2 and so m must be renormalized in the 
second term to compensate the renormalization of X in 
the first term. 

The recursion relations Eqs. (|l9|)-(|2~ij) differ from those 
of Ref. pl| in the m correction to the equation for K. 
This is because of an approximation made there which 
was pointed out to be incorrect by Knops and den Ouden. 
jl6| The approximation, common in the literature, [ ^4|j25| 
is to Taylor expand a cosine term (see the first term in 
Eq. (A8) of Ref. and then to neglect the higher or- 
der terms. As Knops and den Ouden showed, the higher 
order terms are not negligible and in fact contribute in 
a significant way. In our case, the contribution of the 
higher order terms is to give an to correction to the re- 
cursion relation for K that is opposite in sign to that 
obtained by Ichinose and Mukaida. (As discussed in 
the previous paragraph, the sign found here makes more 
physical sense.) However, the m correction to the recur- 
sion relation that most influences the system and the in- 
tegration of the recursion relations turns out to be that 
for y, Eq. (|20|), and this renders the effect of the error 
discussed above less harmful. 

III. ANALYSIS OF THE RECURSION 
RELATIONS 

By integrating Eqs. (]l^)-(|2l|), the analysis of the effects 
of the finite length scale on the 2DCG is possible. One 
sees at once that the size effect, or mass, term is relevant. 
Because y 2 is small, the growth of to with e is nearly 
exponential, which leads to the intrinsic finite size effects 
dominating the integration of the recursion relations at 
large length scales and, in certain temperature regimes, 
at smaller length scales too. In this section we analyze 
this behavior and the temperature regimes where these 
effects are most significant. We will first analyze the RG 
flows, move on to discuss the correlation length, and then 
examine the temperature dependences and the behavior 
of the renormalized stiffness constant. 

In Figure [|, we show a plot the "flows" of the three sys- 
tem parameters as a function of length scale. To facilitate 
the comparison of these results with those of the 2DCG 
in the literature, we use the parameter x — 4/ K — 1. The 
flows are obtained by integrating the recursion relations 
for representative initial values of the three parameters: 
Ki = 8 (Xi = -0.5), m l = 1CT 5 , and Vl = 0.56, 0.58, 
0.60, 0.62, 0.64, 0.66, 0.68. (The initial values of each 
parameter will be denoted below by the subscript i.) The 



recursion relations are integrated until they cease to be 
valid, which is when y or to become large [y(e) > 0.7 
or m(e) > 0.7] and the approximations that we made in 
deriving them break down. As one can see, the flows al- 
ways tend towards large to at the end of the flows. This 
reflects that m (AjJ is a relevant parameter: finite size 
effects will ultimately influence or determine the criti- 
cal behavior. This parameter dominates the end of the 
flows, which corresponds to the largest length scales: fi- 
nite size effects result in free vortices, and the average 
distance between these free vortices corresponds to these 
long lengths. This will become more evident as we ex- 
amine other aspects of the flow. 



m 




FIG. 1. The RG flows for the 2DYG. Because X± is a rel- 
evant parameter, all flows will eventually flow out of the x-y 
plane. The direction of the flows in the figures can be inferred 
from the location of the starting points, as given in the text. 

A projection of these flows onto the x-y plane is shown 
(solid lines) in Figure ^ to better illustrate the effect of m 
on the 2DCG. Also plotted in the figure are the RG flows 
(dashed lines) for the 2DCG (i.e., the to = 2DYG). We 
first explain the RG flows for the 2DCG. The low temper- 
ature (T < Tktb) flows all iterate to y — 0, which corre- 
sponds to zero vortex density since there are no pairs at 
the largest separations for T < Tktb- This follows from 
the relation between vortex fugacity and vortex density. 
The high temperature (T > Tktb) flows tend toward 
y = initially but ultimately towards y — oo. The initial 
decrease reflects the vortex pair density decrease, while 
the upturn and subsequent rise is a result of the increas- 
ing number of free vortices above the transition temper- 
ature. The isotherm that divides the low temperature 
flows from the high temperature flows represents the crit- 
ical temperature. For the linearized recursion relations 
for the 2DCG, the critical isotherm is the line y = —x. 
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FIG. 2. The projection of the RG flows for the 2DYG (solid 
lines) onto the x-y plane. For comparison, the RG flows of 
the 2DCG (dashed lines) are also included. 

Looking now at the RG flows of the 2DYG contrasted 
to the 2DCG, one sees that it is at the largest values of 
e that the largest deviation from the 2DCG flows occurs. 
What is also visible here and in Figure [l] is that none 
of the flows end at y = but all eventually will tend to- 
wards m, y = oo. In other words, all of the flows have the 
same qualitative behavior. This confirms the expected 
behavior that there is no phase transition in the 2DYG 
since there are free vortices present at all temperatures. 
Further, the fact that more RG flows move towards the 
higher temperature limit than for the 2DCG shows that 
the apparent phase transition temperature does shift to- 
ward lower temperature before being wiped out. Along 
with the tendency of the 2DYG flows to move toward 
y = oo, we note that the corresponding tendency for the 
flows to more rapidly approach large x than those of the 
2DCG, which is especially evident for the three highest 
temperatures in Figure ||. This reflects the tendency of 
the renormalized stiffness constant K to approach zero 
more rapidly in the 2DYG. 

It is instructive to look solely at y as a function of 
length scale because this quantity represents the density 
of vortices as a function of length scale. For the same 
temperatures as in Figures l] and |[ we plot y as a func- 
tion of e = hiT in Figure p (solid lines) along with y{e) 
for the 2DCG (dashed lines). For temperatures far below 
Tktb, one can see that the density of vortices is drop- 
ping rapidly with length scale. As mentioned above this 
range of e represents vortex pairs. At large enough length 
scales, there is always an upturn in y(e). This is due 
to free vortices present because of finite size effects. As 
the temperature is increased, while still remaining below 
Tktb, the upturn in y at large e moves to lower values of 
e. At temperatures above Tktb, the behavior is similar 
but more pronounced. This is because the upturn in y 
at large e is always present in the 2DCG for T > Tktb 
since free vortices can be spontaneously created in this 
temperature range. 



FIG. 3. The vortex fugacity y as a function of e, the loga- 
rithm of the length scale for the 2DCG (dashed lines) and the 
2DYG (solid lines) with m, = 10~ 5 . The initial fall is due to 
the decreasing number of vortex pairs as a function of increas- 
ing size. The increase in y at larger length scales represents 
free vortices, which are always present in the 2DYG. 

Further insights can be gained through an examination 
of the correlation length, which is obtained as follows. As 
one integrates the recursion relations, y, K, and to are 
parameterized by e, the natural log of the length scale. As 
mentioned above, the recursion relations are integrated 
until they break down. The value of e at that point is 
denoted by e max . Because there is only one length scale 
in the system and it is the vortex correlation length £(T), 
one makes the association e max — ln£(T). P, p6| , p7| The 
temperature dependence enters through the first integral 
of the recursion relations of the (non-linearized) 2DCG: 
y\ — 2xf + ln(l + Cj) = c where c oc T — Tktb- For small 
to, we assume that this association remains valid. 

The meaning of the 2DCG correlation length below 
Tktb, £-(T), deserves mention. Kosterlitz (2j originally 
defined £_ (T) to be infinite because the susceptibility be- 
low the transition temperature is infinite. Ambegaokar 
et ai, |28| on the other hand, defined a finite diverging 
correlation length for T < Tktb, based on the critical 
behavior of the dielectric constant. (Because the two de- 
scriptions have different meanings, they do not contradict 
one another.) Later, Simkin and Kosterlitz |L3) stated 
that £-(T) is the length scale at which K(e) is essen- 
tially at its asymptotic value In other words, £-(T) 
has meaning both above and below Tktb and here, we 
take it to represent the characteristic size of the largest 
vortex pairs. £+(T), on the other hand, is defined by 
where the recursion relations become invalid and can be 
thought of as the length scale at which the influence of 
free vortices is significant. The relationship of £_ to £ + 
can be seen in Figure [| £- is the length scale at which 
the density of paired vortices becomes small, while £+ is 
the length scale at which the density of free vortices is 
large. In the case of the 2DYG, both £ + and £_ can be 
defined at all temperatures. 
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Defining the correlation length as the quantity £± in- 
troduced by Simkin and Kosterlitz ]T^ | is particularly 
important in the case of the 2DYG. Because m is a rel- 
evant RG parameter, the condition for terminating the 
integration of the recursion relations will always be de- 
termined by m becoming too large. This will result in a 
featureless £ for T < Tktb- But computed in accor- 
dance with Ref. by using the condition that y reaches 
an arbitrary low value, has a clear physical meaning. 
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FIG. 4. e max (the logarithm of the correlation length) ver- 
sus the integration constant c oc T /Tktb — 1 for mt — 1CP 20 , 
10~ 15 , 10~ 10 , and 10~ 7 . The plateau is where the correlation 
length exceeds Aj_. 

The correlation length, defined as explained above, is 
plotted in Figure ||, for various initial values of m, versus 
c oc T — Tktb- As one can see, for each value of mi, 
the correlation length starts to diverge near Tktb but 
then stops growing and reaches a plateau. £ to the left 
of the plateau is determined by y being small: y = 0.01. 
£ in the temperature range of the plateau is determined 
by the large m cutoff, and £ to the right of the plateau 
is determined by the large y cutoff. The abrupt changes 
in £ are due to the cutoff in the integration shifting one 
condition to another. The cutoff value for the small value 
of y is rather arbitrary but the same qualitative behavior 
is found using other values. Furthermore, as one can see 
from Figure [| y = 0.01 does represent a value above 
which the vast majority of the vortices are paired. 

Two temperature quantities can be studied by exam- 
ining the behavior of £ in Figure ^. The first is the tem- 
perature width of the plateau, Tf s . As m is increased, 
Tf s also increases and its dependence on m is well char- 
acterized by 

r /s ocl/ln 2 m. (22) 

It has been previously found that the shift in T c for a 
layered systems has this inverse logarithmic square de- 
pendence upon the strength of the interlayer coupling, 
^9| , p6|j29| indicating a deeper reason for this dependence. 



The second quantity that can be studied, at least quali- 
tatively, is the transition temperature. Before the diver- 
gence of the correlation length is cut off, one can see that 
the effective transition temperature has shifted down. 
This is as expected since the vortex pairs interact more 
weakly. The temperature shift should have the same de- 
pendence as the quantity in Eq. ( |22| ) but this can not be 
verified here because the critical behavior is rounded off. 

We now examine the renormalization of the stiffness 
constant (superfluid density) K. For the 2DCG, this 
quantity, which is evaluated for e = oc, first decreases 
linearly with increasing temperature and then decreases 
more rapidly near Tktb with the temperature depen- 
dence of the well-known square-root cusp. Precisely at 
Tktb, it jumps to zero from its universal value of 4 (in 
these units). 

In the case of the 2DYG, K(oo) is equal to zero, inde- 
pendent of temperature, because there are free vortices at 
all temperatures. Hence, one must look at K for length 
scales £, as shown in Figure |[ £ is the pertinent length 
scale especially for T Tktb where the cutoff deter- 
mining £ is y(e) small, because it characterizes vortex 
pairs, and it is the vortex pair properties (in the pres- 
ence of a current) that allow one to measure K via the 
I-V exponent. K(^) is plotted in Fig. || for various ini- 
tial values of m. One can see for m,; = 0, that K has 
the expected behavior for the 2DCG including the jump 
from universal value of 4 (for our definition of this pa- 
rameter) to zero. As the value of is increased, the 
jump becomes considerably more rounded. The temper- 
ature width over which the jump is smoothed is given by 
Tf s . At the largest value of m,, one can see that even the 
T <c Tktb values of K are depressed. 
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FIG. 5. The renormalized stiffness constant K(£) as a func- 
tion of temperature for m, = 0, 1CT 20 , 10" 15 , 1CT 10 , and 
10" 5 . 

In this work, we have considered the effect of the in- 
trinsic finite size scale Aj_, as opposed to the externally 
imposed finite size of the system. While their effects can 
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be expected to be qualitatively similar and difficult to dis- 
tinguish experimentally (see Ref. ||(| for effects of a finite 
system), there is an important difference. In both cases 
there are free vortices below the transition temperature. 
But the effect of the length X± on the vortex interaction 
is a little larger than that of the system size. As a re- 
sult, we expect that the transition temperature will not 
be depressed as much in the finite extent system. There 
may be more differences in terms of the dynamics, and 
this deserves further study. 

To close this section we remark on the implications 
of this work for the 2D massive SG Hamiltonian itself. 
The recursion relations found here for the 2D massive 
sine-Gordon differ sl ightl y from those of Ref. |p]J. As 
discussed in Section 
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the sign of the correction to 
the recursion relation for K found here is opposite to 
that found by those authors. |21 1 We nonetheless feel that 
the behavior claimed by those authors would be affected 
only quantitatively because, as noted above, it is the m 
correction to the y recursion relation that determines the 
effect of the mass on the system. 



IV. EXPERIMENTAL RAMIFICATIONS AND 
THE EFFECT OF CURRENTS 

In this paper, we have rigorously verified the well- 
known notion that the intrinsic finite size effect char- 
acterized by the relevant parameter Aj_ should eventu- 
ally dominate any critical behavior. One should therefore 
ask why this effect is not more prevalent in experimental 
studies. To address this question, one must incorporate 
the effect of an applied current since most experimental 
studies of KTB behavior in superconductors involve such 
a current. There one expects to see finite size effects for 
T < Tktb when the probing length of the current r c 
(oc 1/1) exceeds A^- H,12| We know from the work here 
that the renormalized value of Aj_ decreases nearly expo- 
nentially, [see Eq. (^f|),1 and so one would expect finite 
size effects to always be seen. However, through a deeper 
examination of this condition, we will see that this is not 
the case. To do that, we must know the behavior of r c 
under renormalization. 

Sujani et al. Q have derived the recursion relations 
for the 2DCG in the presence of a current. These can 
also be obtained by taking the 2D limit of the recursion 
relations of the layered vortex gas in the presence of a 
current, as derived in Ref. ^7j. (See also the work of 
Cserti on 2D dislocation systems and the effect of 
stress.) The latter results include a correction to the 
recursion relation for K due to the current I not found 
in the former. Otherwise, the results are in agreement. 

For the case of the 2DYG in the presence of a cur- 
rent, the recursion relations would consist of those for 
the 2DCG along with the first order corrections from 



both the finite size variable m calculated here and the 
current calculated in Refs. (|] and Refs. (27j. The rele- 
vant recursion relation here is that for /: dl/de = I. 

Let us now return to the condition for seeing inherent 
finite size effects in the presence of a current: r c ^> A^. 
Even though Aj^ decreases nearly exponentially under 
renormalization, it will never get smaller than r c unless 
as its bare value starts out smaller than the bare value 
of r c , because r c is decreasing exponentially under renor- 
malization. This explains the experimental situation. 

We now briefly address the expected effect of finite size 
effects on the current- voltage relationship for the 2DCG. 
The first order effects have been worked out for the lay- 
ered case. [[33) Taking the 2D limit of Equation (11) in 
Ref. ||, (i.e., l\ -> 0,) one obtains §3 



V = I{a + bI 2a ) 1/2 , 



(23) 



where a is the familiar 2DCG I-V exponent and a and 
b are temperature dependent parameters. As expected, 
jl3|,[L4| the I-V relationship is ohmic at small currents, 
due to free vortices induced by finite size effects, and 
then becomes a power-law due to vortex pairs for larger 
currents. The crossover from ohmic to non-ohmic be- 
havior is determined by the vortex correlation length £_ , 
thereby setting limits on a and b. Eq. ( p3f ) is based on 
approximations assuming weak current. We expect that 
a more rigorous treatment of the I-V relation in the pres- 
ence of finite size effects would include extensions to finite 
current, such as a dynamic pd|Jl2] | or finite size scaling 
approach. 



V. SUMMARY 

In this work, we have derived the recursion relations 
for the two-dimensional Yukawa gas, relevant to vortices 
in 2D superconductors. We have examined the renormal- 
ized stiffness constant as a function of length scale, the 
RG flows, and the width of the region around Tktb that 
is no longer critical because of the finite size effect. We 
have also examined the condition for observing finite size 
effect in electrical transport measurements on supercon- 
ducting films. 
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